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Density-functional calculations are used to study various plausible structures of the wurtzite InN 
(0001) and (0001) surfaces. These structures include the unreconstructed surfaces, surfaces with 
monolayers of In or N, several possible coverages and locations of In or N adatoms and vacancies. The 
stable structure of the (0001) surface under N-rich conditions is the unreconstructed, In-terminated, 
surface, while under In-rich conditions the stable surface has a 3/4 monolayer of In atoms. The 
stable structure of the InN (0001) surface corresponds to a full monolayer of In atoms in the atop 
sites (directly above the N atoms) over the entire range of accessible In (or N) chemical potential. 
The atomic structures of the low-energy structures are also discussed. 

PACS numbers: 68.47.Fg 

Keywords: InN surfaces, surface reconstruction, group-Ill nitrides 



I. INTRODUCTION 

Group-Ill nitrides are of considerable interest for ap- 
plications in laser and light-emitting diode n^i^i"^ . Alloy- 
ing different group-Ill elements in nitride compounds 
(e.g., AUGai-ajN or Ina;Gai_a:N) provides a means of 
tuning the band gap over a wide spectral range (e.g., 
from blue to infrared in the case of In^jGai-^^N). The 
growth of high quality group-Ill nitride compounds can 
be a considerable challenge and a barrier to the success- 
ful transition of these materials into applications. The 
connection between vapor phase growth processes and 
film quality can usually be traced back to the interplay 
of growth conditions and surface structure. While there 
have been several theoretical studies on the surfaces of 
group-Ill nitrides, these have largely focused upon AIN 
GaN^i'Pi^iTi^i^iiPi". Relatively little theoretical anal- 
ysis has been devoted to the surfaces of InN, largely be- 
cause the level of interest in this material has been limited 
by difficulty in producing this material. 

InN is an attractive material for electronic and opto- 
electronic applications because it possesses a series of out- 
standing propertiesi^, such as a small effective massif 
(leading to a large electron mobility and a high car- 
rier saturation velocity). The steady-state drift veloc- 
ity in InN is 4.2 x 10^ cm/s is significantly larger than 
that of many III-V competitors, GaN (2.9 x 10'' cm/s), 
AIN (1.7 X 10^ cm/s), and GaAs (1.6 x 10^ cm/s)lij4. 
Since InN exhibits a high peak overshoot velocitjii^, it 
is a good candidate for high-performance heterojunction 
field-effect transistors. 

The growth of high quality InN is difficult compared 
to other Ill-nitrides, such as AIN and GaN, because it 
has a relatively low dissociation temperature and a high 
equilibrium N2 vapor prcssureiS*i^. The lack of suitable 
substrate materials also hinders the use of InN as a het- 
eroepitaxial thin film. Nonetheless, several groupsi2*i2iiS 
have successfully grown good quality InN single crystals. 
Interestingly, they have reported bandgaps of approxi- 
mately 0.9 eV, which is quite small compared to the com- 



monly accepted values of 1.8 cV to 2.1 cV. 

Since InN does not possess a center of symmetry, the 
(0001) and (0001) surfaces are different. The first is 
In-terminated, while the latter is N-terminated. Ex- 
perimental evidence suggests that polarity (termina- 
tion) is an important consideration in the growth of 
high quality group-Ill nitride semiconductors^^. For ex- 
ample, surface with opposite polarity exhibit different 
morphologicsSi. Radio frequency molecular beam epi- 
taxial (rf-MBE) growth of InN on sapphire^S produces 
different polarity surfaces depending on growth temper- 
ature (N-terminated at low temperature, In-terminated 
at high temperature). This demonstrates that growth 
conditions can modify the relative growth rates of sur- 
faces with opposite polarity. The observed change in 
the relative growth rates of surfaces with opposite po- 
larity with temperature may be associated with differ- 
ent activation energies for growth or may be indicative 
a change in surface structure with growth conditions. In 
this work, we determine the equilibrium structure of the 
important (0001) and (0001) surfaces of InN as a function 
of growth conditions. In particular, we examine the role 
of In (or N) chemical potential on the equilibrium sur- 
face structure. We vary chemical potential rather than 
temperature, since the range of chemical potentials that 
are experimentally accessible is much larger than that 
for temperature in the most common group-III nitride 
growth methods (i.e., MOCVD). 



II. THEORETICAL MODELS 

The present calculations are performed within the lo- 
cal density approximation (LDA) framework, as imple- 
mented within the Vienna Simulation Package (VASP)^, 
employing ultrasoft pseudopotentials^^. The important^i 
4d states for In are included as valence electrons in the 
pseudopotentials. Throughout this work, we employed a 
relatively high cutoff energy of 435 eV. We first optimize 
the bulk wurtzite InN structure, yielding the following 
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lattice parameters'^: a = 3.518 A, c = 5.690 A, and 
the internal parameter u = 0.379. These lattice param- 
eters arc in very good agreement with those found from 
experimentii; a = 3.5365 A and c = 5.7039 A. 

For the surface calculations, our supercells contained 
six (0001) bilayers in which the lower four bilayers were 
fixed in the bulk configuration. The upper two bilayers 
and any adatoms or adlaycrs were allowed to relax. Our 
system is larger than those used by many others in stud- 
ies of other group-III nitride surfaces, where typically 
four bilayersSiLSiiS were employed. To prevent unphysi- 
cal charge transfer between the top and bottom slab sur- 
faces, pseudo-hydrogens with fractional charges are used. 
For the In-terminated (0001) surface, the dangling bonds 
in the N layer are saturated with pseudo-hydrogens with 
a fractional charge of 3/4. For the N-tcrminatcd (0001) 
surface, we saturate the dangling bonds in the In layer 
with pseudo-hydrogens with a fractional charge of 5/4. 
We employ a 2 x 2 x 1 Monkhost-Pack set sampling 
scheme for the surface supercell, which gives two irre- 
ducible fc-points. The structural optimizations were ter- 
minated when the magnitude of the HcUmann-Feynman 
force on each ion is less than 1 mRy/Bohr (25.7 meV/A). 
We perform calculations with vacuum thicknesses of both 
12.9 A and 15.7 A. These two sets of calculations yield 
almost identical results, showing the appropriateness of 
the vacuum thickness employed. Hereafter we present 
the results obtained with a vacuum thickness of 15.7 A. 

We consider various vacancy, adatom, trimer, and 
monolayer surface structures that satisfy the electron 
counting rule^. It is useful to refer to Fig. ^a), which 
shows two (0001) layers of atoms in a 2 x 2 surface unit 
cell, in order to understand the structures. The sur- 
faces considered include a perfect surface [as shown by 
the black circles in Fig. ^a)], a monolayer of In (or N) 
atoms directly above each In (or N) terminating atom 
[i.e., this is 4 atop atoms in the 2x2 cell of Fig.^a)], a 
quartcr-monolaycr of adatoms with the adatoms directly 
atop one of the In (or N) atoms (T4 sites) in the 2x2 
cell, a quarter-monolaycr of adatoms with the adatoms 
centered between three surface atoms (113 sites), a single 
vacancy in the terminating plane of the 2x2 cell, and a 
triangle of three In (or N) atoms on H3 or T4 sites (i.e., 
a trimer). The relaxed structure of a representative case 
of an In trimer on the T4 sites is shown in Figs.^b) and 

Me). 

The relative formation energies of the different sur- 
faces are calculated within the thcrmodynamically al- 
lowed range of the In chemical potential ^in. The chem- 
ical potentials of In and N are related through2& 

Wn + A^N = MlnN(bulk)- (1) 

Operationally, we can define the chemical potential 
A*inN(buik) as the total energy per cation-anion pair 
£'(InNbuik) in wurtzite InN at K. The chemical poten- 
tial of In is bounded above by that of the pure, tetragonal 
In metalSi, denoted by = E{In). Similarly, the chem- 
ical potential of N is bounded above by that of a nitrogen 




FIG. 1: (Color online) (a) A schematic of the atom positions 
in the top two atomic planes of a 2 x 2 unit cell of the {0001} 
surface. The larger and smaller filled circles denote indium 
atoms on the surface and nitrogen atoms one plane down, 
respectively, for the In-terminated (0001) surface. The same 
schematic describes the (0001) surface, but with the elemen- 
tal identity of the atoms switched. The locations of the H3 
(hollow) and T4 (top) sites are labeled by dotted circles, (b) 
A view of an In trimer residing on the T4 sites from 
above. The gray and black atoms denote the In and 
N atoms, respectively. Notice that a nitrogen atom 
is exposed in the lower right corner of the image, (c) 
The side view of the 2x2 unit cell. Four pseudo- 
hydrogen atoms (shown in white) are found at the 
bottom of the cell. 



molecule N2, denoted by /.(^ = ^E(N2)- The enthalpy of 
formation of InN is simply 

AH = ^lnN(bulk) ~ A^N ~ Min (2) 

where E{ln), E(N2), and i?(InNbuik) are determined 
from total energy calculations within LDA. The calcu- 
lated enthalpy of formation of InN, AH = —0.376 eV, 
agrees reasonably well with the experimental value22i2S 
of —0.28 eV. The fact that these calculations suggest that 
InN is more strongly bound than in experiments is not 
surprising since the LDA tends to lead to ovcrbinding. 
The small enthalpy of formation implies that InN has a 
low thermal stability compared with other group-III ni- 
trides such as GaN and AIN where the experimental for- 
mation enthalpies22i2S are —1.23 eV (i.e., 4 times larger 
than for InN) and —3.21 eV (i.e., 11 times larger than 
for InN), respectively. 

The formation energy of any structure is calculated 
from 

= Eopt ~ fJ-Ininm - «n) - MlnN(bulk)'^N, (3) 
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FIG. 2: The relative formation energies for the 2x2 unit cell 
of Fig. ^a) as a function of the In chemical potential for the 
(a) (0001) and (b) (0001) surfaces. The formation energies 
are all referenced to that of the surface with a single vacancy 
on the surface atomic plane. 



where i?opt is the total energy of the structure result- 
ing from the optimization of the structure with respect 
to atomic positions and njn and riN arc the number of 
In and N atoms in the structure, respectively. Wc have 
used Eq. to eliminate fj,^ to arrive at the second line 
of Eq. Q • Equations ^ and Q show that the thermo- 
dynamically allowed range of the chemical potential of In 
is Ai/ + /i?„<W„<M?„. 



III. THE (0001) SURFACE 

The relative formation energies for the 2x2 InN (0001) 
unit cell of Fig.^a) are shown in Fig.[2Ia) for the surface 
structures described above as a function of the In chem- 
ical potential. In this figure, we reference the energy to 
that of the 2 X 2 InN (0001) imit cell with a single In 
vacancy in the surface atomic plane. This simplifies the 
comparison of energies between the different structures, 
as well as leads to the cancellation of the contribution to 
-Eopt in Eq. Q from the chemical potential of pseudo- 
hydrogen atoms. We omit the curves for the structures 
of energies greater than 0.25 eV. Note that the slopes of 
the lines in Fig. [2| simply follow from the fact that the 
different structures have different numbers of In and N 
atoms (see Eq. ©). 



The surface with an N adatom on the H3 site has a 
lower energy than that with an N adatom on the T4 
site by 0.70 eV/2 x 2. This is largely due to a stronger 
electrostatic repulsion between the T4 nitrogen adatom 
and the nitrogen atom below. Similar behavior has been 
observed for the AIN^ (0001) and GaNi (0001) surfaces. 
The surface with an In adatom on the T4 site has a lower 
energy than that with an In adatom on the 113 site since 
there is an electrostatic attraction between the In adatom 
on the T4 site and the N atom below. The surface with 
a full In monolayer with atoms on the atop sites has a 
relatively high energy compared with the other (0001) 
surfaces. 

Under N-rich conditions (low In chemical potential), 
the bare 1x1 (In terminated) surface [as shown in 
Fig. n^a)] has the lowest energy. Our calculations also 
show that surfaces in which the In trimer resides on ei- 
ther the T4 or 113 sites yield essentially the same en- 
ergy, despite the fact that In atoms are directly above 
N atoms in the T4 In trimer model. Similar behavior is 
observed for the AlN(OOOl)^ surface where Al trimcrs on 
the T4 or 113 sites have the same energy. This shows that 
arguments based solely on electrostatic interactions are 
not sufficient to explain even the qualitative behavior of 
many of the adatom reconstructed surfaces. Figure |21a) 
shows that under In-rich (large In chemical potential) 
conditions, the stable (0001) surface have In trimers, lo- 
cated on either T4 or H3 sites. The present results for 
the stable (0001) InN surface arc significantly different 
than the same surface in other group-Ill nitrides across 
the entire range of chemical potentials. For AlN(OOOl)^ 
and GaN(0001)i surfaces, the surface with an N adatom 
on the 113 site is found to be energetically favorable un- 
der N-rich conditions, while the surface structure with a 
metal atom (i.e., Al or Ga) on the T4 site is energetically 
favorable under metal-rich conditions. 

Examination of the unreconstructed InN (0001) sur- 
face [as shown in Fig. da)] shows that the In-N bilayer 
spacings for the first and second bilayers are 0.654 A and 
0.655 A, respectively. These are slightly smaller than the 
ideal, bulk spacing of 0.688 A. The spacing between In- In 
atomic planes are 2.862 A and 2.825 A for the ffist and 
second pairs, as compared with a bulk spacing of 2.845 A. 
This shows that the spacing between the first In-In layer 
expands, while the second contracts, relative to the bulk 
spacing. Similar oscillatory surface relaxations are ob- 
served in other systems (e.g., see Ref. ^3)- 

For the InN (0001) surface with a single In vacancy on 
the surface layer, wc find that the top-layer In atoms un- 
dergoes sp^ bonding configuration by relaxing vertically. 
This reduces the bond length between the threefold- 
coordinated In and the threefold-coordinated N atoms 
from the ideal value of 2.145 A to 2.050 A. This 4.4 % 
reduction in length is comparable to that for the case of 
AIN^, where a 6% contraction in the distance between 
the Al and N atoms in the vacancy model of AlN(OOOl) 
surface was observed. The vertical movement of the top 
layer atoms reduces the top In-N bilayer spacing from the 
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ideal value of 0.688 A to 0.350 A. 

The surface for with an In trimer on the T4 or H3 sites 
have the same energy and are the most stable surfaces 
at large In chemical potential, as discussed above. It is 
interesting to note that the interlayer spacing between 
the trimer In atoms and the In atoms in the next plane 
are 2.61 A for the T4 site trimer and 2.60 A for the H3 
site trimer. Hence, these two surface configurations are 
very similar both structurally and energetically. 



IV. THE (0001) SURFACE 

The relative formation energies for the 2x2 InN (0001) 
unit cell of Fig.QIa) are shown in Fig.|3b) for the surface 
structures described above as a function of the In chem- 
ical potential. In this figure, we reference the energy to 
that of the 2x2 InN (OOOT) unit cell with a single N 
vacancy in the surface atomic plane. No results are pre- 
sented for the surfaces with N adatoms or N trimers since 
these systems are unstable. Northrup et al^ argued that 
the relatively small size of the N atom makes surfaces 
with N adatoms or trimers on the AlN(OOOl) highly un- 
stable. It is therefore not surprising that our calculations 
also show that N trimer or N adatom models are not sta- 
ble since the single-bond covalcnt radius^ of In (1.48 A) 
is considerable larger than that of Al (1.30 A). 

The relaxed unreconstructed InN (0001) surface (N- 
terminated) has a very large formation energy, relative 
to most of the other surface structures. We find that 
the surfaces with an In adatom on the 113 site has an 
energy that is 1.83 eV/(2 x 2) lower than that with the 
In adatom on the T4 site. This may be attributed to 
the fact that there is greater electrostatic repulsion be- 
tween the T4 In adatom and the In atom directly below 
it than between the 113 In adatom and the more distant 
In atoms two atomic planes below. Similar argument also 
holds for In trimer models, where the H3 In trimer has a 
lower in energy than the T4 In trimer by 1.23 eV/(2 x 2). 
The surface with an In monolayer with the atoms on the 
atop sites has the lowest energy across the entire allowed 
range of In chemical potential. The next lowest energy 
structure under N-rich (low In chemical potential) con- 
ditions is surface with a single In adatom in the 2x2 
unit cell on the H3 site. These results contrast with the 
results for other group-Ill nitrides (0001) surfaces. In 
AIN and GaN, the lowest energy surface structure corre- 
sponds to an Al monolayer (atop sites occupied)^ or Ga 
monolayer™ only under metal-rich conditions. However, 
under N-rich conditions, the (0001) surfaces with a single 
Al or Ga atom on the H3 site in the 2x2 unit cell have 
the lowest energy. 

Examination of the InN (0001) surface with an In 
monolayer with the In atoms on the atop sites is the 
lowest energy structure over the entire chemical poten- 
tial range, as mentioned above. The In-N bilayer spacing 
for the top two bilayers are 0.519 A and 0.510 A, as com- 
pared with a bulk spacing of 0.668 A. The In-In atomic 



plane separations are 2.482 A, 2.423 A, and 2.593 A for 
the first three pairs (counted from the surface), as com- 
pared with the bulk value of 2.845 A. This is a substan- 
tial contraction that oscillates as it decays slowly into the 
bulk. 

The surface structure corresponding to an In adatom in 
the H3 site also has a relatively low energy. This adatom 
is situated 0.808 A above the triangle of three N atoms in 
the plane below. The corresponding In-N bond length is 

I. 980 A. The three N atoms that form this triangle (i.e., 
nearest neighbors of the In H3 adatom) are displaced 
upward, toward the In adatom by 0.424 A relative to the 
other, threefold-coordinated surface N atoms. 

In the InN(OOOl) surface structure with a single N- 
vacancy in the 2x2 unit cell, the top-layer N atoms un- 
dergo substantial vertical relaxation. In this structure, 
the N-In bilayer spacing is 0.328 A as compared with 
the bulk spacing of 0.688 A. The bond length between 
the threefold-coordinated N and In atoms is also reduced 
from the ideal value of 2.145 A to 1.898 A. This is an 

II. 5% reduction in bond length, which is much larger 
than the 6% contraction observed in the N vacancy sur- 
face structure of AlN(OOOl)^. This large contraction may 
be attributed to the system trying to move toward a con- 
figuration that is sp^ bonded. 



V. CONCLUSIONS 

We have performed density-functional calculations to 
determine the stable structure of both the InN (0001) and 
(0001) surfaces as a function of the In (or N) chemical 
potential. We considered twelve plausible surface struc- 
tures. Several earlier studies considered the structure 
of this surface in other group-Ill nitrides. While struc- 
turally similar, InN is quite distinct, as indicated by a 
much lower heat of formation than the others (BN, AIN, 
GaN). 

For the (0001) surface, we find that the relaxed, un- 
reconstructed InN (0001) (In terminated) surface is sta- 
ble under N-rich (low In chemical potential) conditions. 
On the other hand, under In rich (large In chemical po- 
tential) conditions, the stable structure corresponds to 
a trimer of In adatoms on either H3 or T4 sites (this 
is 3/4 of an In monolayer). These stable (0001) surface 
structures differ from those found for two other closely 
related group-Ill nitrides over the entire range of accessi- 
ble chemical potentials. The stable structure of the InN 
(0001) surface corresponds to a monolayer of In atoms 
directly above the surface N atoms. This structure is 
stable over the entire range of accessible chemical poten- 
tials. In AIN and GaN, the same surface structure was 
only found to be stable under metal-rich conditions. 

The present results provide guidance for the exploita- 
tion of different thin film growth methods / conditions in 
order to exploit the dependence of equilibrium recon- 
structions on chemical potential. Unlike in GaN or AIN, 
the stable structure of the important (0001) surface can 
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be manipulated by changes in growth condition. 



supported by Visiting Investigator Program, Agency for 
Science, Technology and Research (A*STAR), Singapore. 



VI. ACKNOWLEDGMENTS 



The authors gratefully acknowledge useful discussions 
with Y. Y. Sun on surface calculations. This work was 



^ S. Nakamura, M. Senoh, S. Nagahama, N. Iwasa, T. Ya- 
mada, T. Matsushita, H. Kiyoku, and Y. Sugimoto, Jpn. 
J. Appl. Phys. 35, Part II, L74 (1996). 
S. Nakamura and G. Fasol, The Blue Laser Diode: GaN 
Based Light Emitters and Lasers (Springer, Berlin, Heidel- 
berg, 1997). 

^ Y. Taniyasu, M. Kasu, and T. Makimoto, Nature 441, 325 
(2006). 

* J. E. Northrup and J. Neugebauer, Phys. Rev. B 53, 
R10477 (1996). 

^ J. E. Northrup, R. DiFelice, and J. Neugebauer, Phys. Rev. 

B 55, 13878 (1997). 
" A. R. Smith, R. M. Feenstra, D. W. Greve, J. Neugebauer, 

and J. E. Northrup, Phys. Rev. Lett. 79, 3934 (1997). 

K. Rapcewicz, M. BuongiornoNardelli, and J. Berholc, 

Phys. Rev. B 56, R12725 (1997). 
® J. Fritsch, O. F. Sankey, K. E. Schmidt, and J. B. Page, 

Phys. Rev. B 57, 15360 (1998). 
^ F. H. Wang, P. Kriiger, and J. Pollmann, Phys. Rev. B 64, 

035305 (2001). 

^° C. D. Lee, Y. Dong, R. M. Feenstra, J. E. Northrup, and 

J. Neugebauer, Phys. Rev. B 68, 205317 (2003). 
" V. Timon, S. Brand, S. J. Clark, M. C. Gibson, and R. A. 

Abram, Phys. Rev. B 72, 035327 (2005). 

A. G. Bhuiyan, A. Hashimoto, and A. Yamamoto, J. Appl. 

Phys. 94, 2779 (2003). 

S. N. Mohammad and H. Morkoc, Prog. Quant. Electr. 20, 
361 (1996). 

" B. E. Foutz, S. K. O'Leary, M. S. Shur, and L. F. Eastman, 
J. Appl. Phys. 85, 7727 (1999). 

J. B. MacChesney, P. M. Bridenbaugh, and P. B. 
O'Connor, Mater. Res. BuU. 5, 783 (1970). 
O. Ambacher, M. S. Brandt, R. Dimitrov, T. Metzger, 
M. Stutzmann, R. A. Fischer, A. Miehr, A. Bergmaier, 
and G. DoUinger, J. Vac. Sci. Technol. B 14, 3532 (1996). 
" V. Y. Davydov, A. A. Klochikhin, R. P. Seisyan, V. V. 



Emptsev, S. V. Ivanov, F. Bechstedt, J. Furthmiiller, 
H. Harima, A. V. Mudryi, J. Aderhold, O. Semchinova, 
and J. Graul, Phys. Stat. Sol. (B) 229, Rl (2002). 
J. Wu, W. Walukiewicz, K. M. Yu, J. W. A. HI, E. E. 
Haller, H. Lu, W. J. Schaff, Y. Saito, and Y. Nanishi, Appl. 
Phys. Lett. 80, 3967 (2002). 

Y. Nanishi, Y. Saito, and T. Yamaguchi, Jpn. J. Appl. 
Phys. 42, Part I, 2549 (2003). 

Y. Saito, Y. Tanabe, T. Yamaguchi, N. Teraguchi, 

A. Suzuki, T. Araki, and Y. Nanishi, Phys. Stat. Sol. (B) 
228, 13 (2001). 

B. Daudin, J. L. Rouviere, and M. Arlery, Appl. Phys. 
Lett. 69, 2480 (1996). 

G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

2^ D. Vanderbilt, Phys. Rev. B 41, R7892 (1990). 
^'^ A. F. Wright and J. S. Nelson, Phys. Rev. B 51, 7866 
(1995). 

^■"^ C. K. Can, Y. P. Feng, and D. J. Srolovitz, Phys. Rev. B 
73, 235214 (2006). 

G. X. Qian, R. M. Martin, and D. J. Chadi, Phys. Rev. B 
38, 7649 (1988). 

J. Graham, A. Moore, and G. V. Raynor, J. Inst. Metals 
84, 86 (1956). 

W. A. Harrison, Electronic structure and the properties of 
solids (Dover Publications, New York, 1989). 

C. Kittel, Introduction to Solid State Physics (John Wiley 
and Sons, New York, 1996), 7th ed. 

Y. Y. Sun, H. Xu, J. C. Zheng, J. Y. Zhou, Y. P. Feng, 
A. C. H. Huan, and A. T. S. Wee, Phys. Rev. B 68, 115420 
(2003). 

F. A. Cotton, G. Wilkinson, and P. L. Cans, Basic Inor- 
ganic Chemistry (John Wiley and Sons, New York, 1995) , 
3rd ed. 

C. K. Can and D. J. Srolovitz (2006), unpublished. 



